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G q -coupled plasma membrane receptors activate phospholipase C (PLC), which hydrolyzes membrane phospha- 
tidylinositol 4,5-bisphosphate (PIP2) into the second messengers inositol 1,4,5-trisphosphate (IP3) and diacyl- 
glycerol (DAG). This leads to calcium release, protein kinase C (PKC) activation, and sometimes PIP 2 depletion. 
To understand mechanisms governing these diverging signals and to determine which of these signals is respon- 
sible for the inhibition of KCNQ2/3 (K v 7.2/7.3) potassium channels, we monitored levels of PIP2, IP3, and calcium 
in single living cells. DAG and PKC are monitored in our companion paper (Falkenburger et al. 2013. /. Gen. 
Physiol, http://dx.doi.org/10.1085/jgp.201210887). The results extend our previous kinetic model of G q -coupled 
receptor signaling to IP 3 and calcium. We find that activation of low-abundance endogenous P2Y 2 receptors by a 
saturating concentration of uridine 5 '-triphosphate (UTP; 100 uM) leads to calcium release but not to PIP 2 deple- 
tion. Activation of overexpressed Mj muscarinic receptors by 10 uM Oxo-M leads to a similar calcium release but 
also depletes PIP 2 . KCNQ2/3 channels are inhibited by Oxo-M (by 85%) , but not by UTP (<1 %) . These differences 
can be attributed purely to differences in receptor abundance. Full amplitude calcium responses can be elicited 
even after PIP 2 was partially depleted by overexpressed inducible phosphatidylinositol 5-phosphatases, suggesting 
that very low amounts of IP 3 suffice to elicit a full calcium release. Hence, weak PLC activation can elicit robust 
calcium signals without net PIP 2 depletion or KCNQ2/3 channel inhibition. 
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INTRODUCTION 

This and our companion paper in this issue (Falkenburger 
et al.) concern quantitative properties of signaling by 
seven-transmembrane receptors of the plasma mem- 
brane coupling to G q ("G q PCR") , which are examined 
experimentally and interpreted by kinetic modeling. 
Originally motivating this work was our long-standing 
interest in the regulation of the KCNQ2/3 potassium 
ion channel through activation of G q PCRs. These recep- 
tors mediate diverse, important responses of neuronal 
and nonneuronal cells to stimuli such as light sensation 
in Drosophila melanogasler photoreceptors, contraction 
of vascular myocytes in response to adrenaline, behav- 
ioral and mood regulation in response to serotonin 
and dopamine, the effects of cholinergic drugs treating 
Alzheimer's disease, and the inhibition of KCNQ2/3 
(K v 7.2/7.3) potassium channels. 

G q activates PLC, which cleaves the plasma membrane 
phospholipid phosphatidylinositol 4,5-bisphosphate 
(PIP 2 ) into the cytosolic messenger inositol 1,4,5-trisphos- 
phate (IP3) and membrane-bound diacylglycerol (DAG) . 
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IPs binding to IP3 receptors (IP 3 Rs) at the ER triggers 
the release of calcium into the cytosol. Calcium and DAG 
activate PKC. In addition, PLC may deplete its substrate, 
PIP 2 , which is an activating cofactor for many ion chan- 
nels and other membrane proteins (Suh and Hille, 
2008; Logothetis et al., 2010). We explore all these ac- 
tions here. 

Intriguingly, the consequences of PIP 2 hydrolysis — 
(a) calcium release, (b) PKC activation, (c) net PIP 2 de- 
pletion, and (d) channel inhibition — do not always occur 
together. Moreover, most cells harbor several different 
kinds of receptors coupled to G q , whose signaling re- 
sponses sometimes differ. For example, in sympathetic 
neurons, bradykinin (through B 2 receptors) and ace- 
tylcholine (through Mi muscarinic receptors [MiRs]) 
stimulate production of DAG, but only bradykinin leads 
to calcium release (Delmas and Brown, 2002), and only 
acetylcholine leads to observable PIP 2 depletion (Zaika 
et al., 2011). The existence of signaling microdomains 
and the calcium-induced acceleration of PIP 2 synthesis 
have been identified as factors contributing to this dis- 
sociation in sympathetic neurons (Delmas and Brown, 
2002; Zaika etal., 2011). 
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We believe that knowing the quantitative requirements 
for PKC activation, calcium response, and PIP 2 deple- 
tion will be valuable for understanding how such signal- 
ing specificity can work, how PLC activation can lead to 
calcium release but not PIP 2 depletion, and how the 
effects of PLC activation through IVLR can differ from 
those of PLC activation through other G q PCRs. There- 
fore, we performed a quantitative analysis of signaling 
events downstream of PLC, using fluorescent reporters 
for PIP2, IP3, and calcium, and in Falkenburger et al. 
(2013) for DAG and PKC. In that paper, these measure- 
ments are then used to extend a kinetic model of PLC 
signaling that builds on earlier versions (Horowitz et al., 
2005; Falkenburger et al., 2010a,b). We show, for in- 
stance, that the requirements for calcium release and 
KCNQ2/3 channel inhibition are quite different: cal- 
cium release can be evoked with a far lower density of 
activated receptors. 

MATERIALS AND METHODS 

Cell culture and plasmids 

tsA-201 cells were cultured in DMEM (Gibco) with 10% serum and 
0.2% penicillin/ streptomycin and passaged every 5 d. Cells were 
transiently transfected at ^75% confluency with Lipofectamine 
2000 (10 pi for a 3-cm dish; Invitrogen) and 0.5-1.2 pg DNA per 
plasmid. Cells were transfected 1 or 2 d before photometry experi- 
ments and 2 d before patch-clamp experiments. Cells were plated 
on polylysine-coated glass chips 12 h before experimentation. 

The following plasmids were used: dark (nonfluorescent) 
and eYFP-labeled mouse MjR (IV^R and MiR-YFP; provided by 
N. Nathanson, University of Washington, Seattle, WA); human 
purinergic receptor P2Y 2 (P2Y 2 R; The Missouri S&T cDNA Re- 
source Center); EPAC1 (exchange protein directly activated by 
cAMP) Forster resonance energy transfer (FRET) probe for cAMP 
(provided by M. Lohse, University of Wiirzburg, Wurzburg, 
Germany) ; human KCNQ2 (provided by D. McKinnon, State Uni- 
versity of New York, Stony Brook, NY) ; human KCNQ3 (provided by 
T. Jentsch, Leibniz-Institut fur Molekulare Pharmakologie, Berlin, 
Germany); eCFP-PH(PLCSl), eYFP-PH(PLCBl), and CFP-CAAX 
(from K-Ras; provided by K. Jalink, The Netherlands Cancer Insti- 
tute, Amsterdam, Netherlands); the zebrafish voltage-sensitive 
phosphatase Dr-VSP-IRES-GFP (Dr-VSP; provided by Y. Okamura, 
Osaka University, Osaka, Japan) ; "Dark" Dr-VSP (Falkenburger et al., 
2010b); and the IP 3 reporters IRIS-1 (provided by K. Mikoshiba, 
Institute of Physical and Chemical Research Brain Science Insti- 
tute, Wako, Japan) and LIBRA version III (LIBRAvIII; provided 
by A. Tanimura, Health Sciences University of Hokkaido, Tobetsu, 
Japan) . IRIS-1 is cytosolic, whereas LIBRAvIII is membrane local- 
ized by palmitoylation (membrane-targeting sequence of GAP43). 
Version III of LIBRA has the pH-stable YFP mutant Venus instead 
of eYFP and derives from rat IP3R type III. The small interfering 
RNA (siRNA) for P2Y 2 was purchased from Santa Cruz Biotech- 
nology, Inc. Membrane-targeted FRB (LDR) and CFP-FKBP-Inp54p 
were generated by T. Inoue (Johns Hopkins University, Baltimore, 
MD) and described previously (Suh et al., 2006). Hereafter, we refer 
to fluorophores simply as "CFP" or "YFP" regardless of whether 
regular or enhanced fluorescent proteins were used. 

Electrophysiology 

KCNQ2/3 currents were recorded in whole-cell gigaseal voltage- 
clamp or perforated-patch configuration using borosilicate glass 
pipettes with a resistance around 2.2 MCI. Internal solution was 



(mM): 175 KC1, 5 MgCl 2 , 5 HEPES, 0.1 K4BAPTA, 3 Na 2 ATP, and 
0.1 Na 3 GTP, pH 7.4 (KOH). Recordings used an EPC9 amplifier 
with Patchmaster 2.35 software (HEKA) . Currents were filtered at 
2.9 kHz. Sample intervals were 200 ps or slower. Series resistance 
was compensated by 70% after compensation of fast and slow 
capacitance. Leak was not subtracted. Holding potential was 
— 60 mV. KCNQ2/3 current was quantified by measuring tail cur- 
rents. Every 2 s, the membrane was depolarized to —20 mV for 
400 ms and repolarized to —60 mV. KCNQ2/3 current activates 
slowly upon depolarization and deactivates slowly upon repolar- 
ization (see Fig. 1 A, bottom right). KCNQ2/3 tail currents were 
measured by comparing current at 20 and 400 ms after repolar- 
ization to — 60 mV. 

For perforated-patch recordings, amphotericin B was dissolved 
in DMSO on the day of recording (1 mg/10 pi), sonicated, diluted 
in internal solution to a working concentration of 300-500 pg/ ml, 
and sonicated again. The tip of the patch pipette was filled 
with amphotericin B-free Ringer's solution by capillary action. 
After <5 min, access resistance was usually low enough to record 
KCNQ2/3 current. 

Perfusion, extracellular buffer, and temperature 

Cells were recorded in a 100-pl chamber continuously superfused 
(1 ml/min) with Ringer's solution containing (mM): 160 NaCl, 
2.5 KC1, 2 CaCl 2 , 1 MgCl 2 , 10 HEPES, and 8 glucose, pH 7.4 
(NaOH). All drugs were applied in the superfusate. As a mea- 
sure of exchange time, when 2.5 mM KC1 Ringer's solution was 
replaced by 30 mM KC1 Ringer's solution (shifting the K + rever- 
sal potential positive of —60 mV), KCNQ2/3 tail currents were 
inverted with a delay of <1 s and a time constant of 2 s, which 
is sufficiently fast for the kinetics investigated here. All mea- 
surements reported in this paper were done at room tempera- 
ture, 21-23°C. 

Photometric calcium measurements 

All optical measurements of calcium and FRET used a monochro- 
matic light source (Polychrome IV; TILL Photonics) for epifluo- 
rescence illumination and one or two TILL photodiodes with 
suitable filters and dichroic mirrors for photometry. Photo- 
diode voltages were acquired by Patchmaster (sampling 200 ps 
or slower) , and after recording from each cell, an area of the cov- 
erslip without cells was measured as background. Cytosolic free 
calcium was measured with the low-affinity indicator Fura-4F 
(Invitrogen). For measurements without electrophysiology or 
in perforated-patch configuration, cell-permeable Fura-4F-AM 
ester was diluted to 2 pM in Ringer's solution supplemented with 
0.2% pluoronic F-68. Cells were loaded for 40 min and then incu- 
bated in regular Ringer's solution for an additional 30 min to 
allow complete de-esterification of AM esters. For measurements 
in whole-cell configuration, 0.1 mM of cell-impermeable Fura-4F 
salt was added to the intracellular solution in the patch pipette. 
Fura-4F fluorescence was measured by stepping the excitation 
light to 340 nm for 100 ms and 380 nm for 20 ms every 4 s. (The 
recording time with 340-nm excitation was longer than for 380-nm 
excitation to compensate for the poorer signal-to-noise ratio with 
340-nm excitation.) Excitation light was reflected by a 415-dclp 
or a three-color (89006bs; Chroma Technology Corp.) dichroic 
mirror. Fluorescence was detected using one photodiode with 
a 535/30-nm emission filter. Background fluorescence was sub- 
tracted, and the ratio of emission with 340-nm excitation (F 340 ) to 
emission with 380-nm excitation (F 380 ) was calculated offline using 
a custom macro for IGOR Pro 6.0 (WaveMetrics) . 

The following features were extracted from traces of Fura-4F 
ratio (F 34 o/F 380 ): baseline, mean ratio of the points preceding 
agonist application; peak, the maximum ratio reached within 20 s 
after agonist application; time to half-maximum, the delay be- 
tween start of agonist application and the first time point where 
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the ratio exceeded 50% of the distance between baseline and 
peak; duration, the delay between time to half-maximum and 
the first time point when the ratio fell below 10% of the distance 
between baseline and peak. (The level of 10% was chosen to in- 
clude a late "hump" or plateau of the Fura-4F signal in response 
to G q PCR activation.) 

Three methods for obtaining Fura-4F calibration curves were 
compared: (1) adding Fura-4F salt to solutions with known cal- 
cium concentrations and measuring droplets of these solutions, 
(2) measuring cells into which such solutions had been dialyzed 
by the whole-cell pipette, and (3) superfusing cells loaded with 
cell-permeable Fura-4F/AM with such solutions in the presence 
of 15 uM ionomycin. These calibrations yielded similar results. The 
calcium concentration corresponding to a given r = F 340 /F 380 is 
given by: [calcium] =K' (r - r min )/(r max - r), where r min is F 34 o/F 380 
with 20 mM EGTA and r max is F 34 o/F 38 o with 2 mM calcium. Values 
of rmin and rmax were dependent on the dichroic mirror used in 
the microscope. The "triple dichroic" we used to measure FRET 
(89006bs; Chroma Technology Corp.) reflects short wavelength 
excitation light slightly less well than a dichroic mirror designed 
specifically for Fura measurements (DCLP 415). Accordingly, rmin 
and rmax were 0.017 and 0.245 for the 89006bs, and 0.050 and 
0.737 for the DCLP 415. The calcium concentration at which the 
F340/F380 ratio change is half-maximal was K' = 30 uM, regard- 
less of the dichroic used. 89006bs was used for experiments as 
in Fig. 7 C, and DCLP 415 was used in experiments as in Figs. 1 A, 
2, 4, and 7 A. 

We chose a low-affinity calcium indicator to better discriminate 
higher levels of cytosolic free calcium. Using such an indicator 
means that it is difficult to resolve low levels of resting calcium prop- 
erly. It is for this reason that the figures show fura ratio F 34 o/F 380 . 

Photometric FRET measurements 

Epifluorescence photometry measured FRET from CFP to YFP 
using a three-color dichroic mirror in the microscope, two pho- 
todiode detectors, and the excitation wavelength scanned in a 
ramp as described in Falkenburger et al. (2010b). Unless otherwise 
noted, sweeps were repeated every 2 s. The raw fluorescence data 
were corrected for background and for bleedthrough of CFP 
emission into the YFP recording channel, yielding the corrected 
fluorescence values CFP C (480/ 40-nm emission with 440-nm exci- 
tation) and YFP C (535/30-nm emission with 440-nm excitation). 
As before, FRET was expressed as the ratio FRETr = YFP C / CFP C . 
CFP C , YFP C , and FRETr are in arbitrary, apparatus-dependent units. 
In two panels as indicated, the points were mildly smoothed with a 
binomial filter that weighted immediate neighbors by a 0.25/ 
0.5/0.25 algorithm. 

For time courses of FRETr responses, we determined the fol- 
lowing: baseline, mean FRETr of the points preceding agonist 
application; amplitude, the difference in FRETr between baseline 
and maximum or minimum for the response, the delay from 
agonist application to the onset of the response, and the expo- 
nential time constant of the rising and falling phases of the signal 
(by the fitting procedure of Jensen et al., 2009) ; and the duration 
at half-maximum. 

Western blot 

Cells were harvested in ice-cold PBS, centrifuged, and resus- 
pended in 50 pi of lysis buffer (PBS with 1% Triton X-100, 1:100 
EDTA, and protease inhibitor cocktail; Thermo Fisher Scientific). 
Lysates were incubated for 30 min on ice, triturating regularly. 
Lysates were cleared by centrifugation (13,000 g for 20 min at 
4°C), and the supernatant was transferred to new tubes. 20 pg 
of protein was separated by SDS-PAGE using standard techniques 
(Falkenburger et al., 2010a) . The primary antibody against human 
P2Y 2 (Santa Cruz Biotechnology, Inc.) was used at 1:500. The sec- 
ondary antibody (1:30,000; goat anti-rabbit; Kirkegaard & Perry 



Laboratories, Inc.) was coupled to horseradish peroxidase and 
visualized by chemiluminescence using a digital imaging system 
(Alphalmager; Alphalnnotech) . Bands were quantified by mea- 
suring the mean intensity in each lane with the same-sized region 
of interest. Results from four independent experiments were 
averaged for display. Results did not change when intensities for 
P2Y 2 were normalized to amounts of actin as a loading control. 

Modeling 

To summarize the results and test mechanistic hypotheses, a ki- 
netic model of G q PCR signaling was formulated and solved in the 
Virtual Cell simulation environment (University of Connecticut 
Health Center) . Results of the model appear in many figures of 
this paper, but the rationale, design, and parameters of the model 
itself are described only in our companion paper (Falkenburger 
etal., 2013). 

Statistics 

Summarized data include one data point per cell. Numbers, bars, 
and markers represent mean ± SEM. The notation n = xx refers to 
the number of cells averaged unless otherwise noted. Compari- 
sons were made by one-way AN OVA and Tukey posthoc tests using 
GraphPad Prism software. 

Online supplemental material 

Fig. SI shows that uridine 5 '-triphosphate (UTP) does not inhibit 
KCNQ2/3 currents in whole-cell configuration. In Fig. S2, en- 
dogenous P2Y receptors change calcium but not cAMP. Fig. S3 
shows expression patterns of FRET probes/pairs. Fig. S4 is a com- 
parison of the IP 3 FRET probes LIBRAvIII versus IRIS-1. Fig. S5 
shows modeling related to IP 3 and LIBRAvIII. In Fig. S6, time 
course of PIP 2 (PH-domain FRETr) during intermittent voltage- 
sensitive 5-phosphatase (VSP) activation is shown. Fig. S7 shows 
calcium rise in response to caffeine and ryanodine. Fig. S8 shows 
the origin of the plateau in the calcium and Fura-4F responses, and 
in Fig. S9, the time course of acceleration of PIP 2 synthesis is 
established. Figs. S1-S9 are available at http://www.jgp.org/cgi/ 
content/full/jgp. 201210886/DC1. 

RESULTS 

Activation of endogenous P2Y receptors evokes calcium 
release but not significant PIP 2 depletion 

To gain insight into the IP3 branch of signaling elicited 
by activation of G q PCR and PLC, we monitored levels of 
PIP9, IP3, and calcium by optical probes and recorded 
KCNQ2/3 potassium currents in single living tsA-201 
cells. The overall design was to compare responses with 
a purinergic agonist to those with a muscarinic agonist 
in cells transfected with MjRs. tsA-201 cells do not ex- 
press endogenous MiRs. We previously measured both 
endogenous and overexpressed levels of components 
of the M 4 R signaling cascade (Falkenburger et al., 
2010a) and estimated the density of overexpressed 
MjRs as being several orders of magnitude higher (500- 
1,000 per pm 2 ) than typical estimates of endogenous 
receptors ( 1 per pm 2 ) such as endogenous purinergic 
receptors. Henceforward, we refer to transfected MjRs 
as "high-density" MiRs. 

First, we measured calcium using Fura-4F. Despite 
the large expected difference in their receptor densities, 
we found that the Fura-4F response was similar when 
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we activated endogenous purinergic receptors using 
a saturating concentration of 100 pM UTP or activated 
high-density MiRs using 10 pM oxotremorine-M (Oxo-M; 
Fig. 1, A and B). In perforated-patch recording with ei- 
ther stimulus, the Fura-4F ratio (F 3 4o/F 38 o) rose to values 
corresponding to 3 pM of free calcium. On average, 
the amplitude of the Fura-4F response to UTP was 98% 
of the Oxo-M response (Fig. 1 C) , and the duration of the 
response (50-10% of maximum) was similar (Fig. 1 D). 
However, the UTP response did rise more slowly, taking 
approximately three times as long to reach half-maximal 
amplitude (Fig. 1 E), suggesting that the activation of 
G q and PLC by UTP was less intense. The ambient rest- 
ing calcium was below the threshold for reliable quanti- 
tation by the low-affinity Fura-4F dye. 

When we recorded KCNQ2/3 potassium currents in 
the perforated-patch configuration simultaneously with 
Fura-4F signals, we were surprised to find that KCNQ2/3 
currents were not inhibited by 100 pM UTP although 
they were inhibited by 10 pM Oxo-M (Fig. 1, A and F). 
The findings were similar in whole-cell configuration 
(Fig. Sl,AandB),with<l% inhibition by UTP (22 cells) 
and 84% inhibition by Oxo-M (17 cells). Under the per- 
forated-patch configuration, peak current was gradually 
growing (spontaneously) in some cells, sometimes lead- 
ing to apparent "negative" values of inhibition by ago- 
nist. We show perforated-patch configuration in Fig. 1 
because UTP did not elicit a calcium response in whole- 
cell configuration (Fig. SI , A and B) . The absence of a 



whole-cell calcium response with UTP could be a conse- 
quence of (a) diffusion of IP 3 and calcium into the patch 
pipette, thereby reducing IP 3 - and calcium-induced cal- 
cium release, and/ or (b) an inhibitory effect of Mg 2+ 
diffusion into the cell (Volpe et al., 1990). 

The contrast between responses of P2Y receptors and 
MiRs thus raises a challenging question: How do two re- 
ceptors produce a full-amplitude calcium signal, whereas 
only one of them inhibits KCNQ2/3 channels? To ad- 
dress this, we first identified the endogenous subtype of 
P2Y receptors in tsA-201 cells. 

P2Y 2 R is the major endogenous G q PCR responsible 
for UTP-activated calcium rises 

tsA-201 cells express P2Y receptor subtypes 1, 2, 11, 12, 
and 14 (Atwood et al., 2011). Of these, only P2Y 2 R and 
P2Y n R are directly activated by UTP. Both use G q as 
their primary transduction pathway, with P2YnR also 
using G s as a secondary transduction mechanism. We 
found that UTP does not produce a detectable change 
in cAMP levels as measured by a FRET sensor based on 
EPAC (Nikolaev et al., 2004; Fig. S2), suggesting that 
P2Y 2 R would more likely be the major endogenous re- 
ceptor responsible for rises in cytosolic calcium. In the 
following experiments, we confirmed the expression 
and functional importance of P2Y 2 Rin tsA-201 cells by 
(a) Western blot analysis, (b) testing the effects of a 
P2Y 2 R antagonist, and (c) P2Y 2 R "knockdown" by an 
established siRNA. 
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Figure 1. Activating endogenous purinergic re- 
ceptors initiates a calcium rise but does not in- 
hibit KCNQ2/3 current. tsA cells transfected with 
MiR and KCNQ channel subunits were loaded 
with Fura4F-AM, fluorescence was recorded by 
photometry, and simultaneously, KCNQ2/3 po- 
tassium currents were recorded in perforated- 
patch configuration. (A; left) Time course of 
Fura-4F ratio and tail-current amplitudes. (Top 
right) Component single sweeps of data collec- 
tion showing photodiode output as wavelength is 
stepped to 340 and 380 nm. (Bottom right) De- 
polarization-evoked KCNQ currents (steps from 
— 60 to —20 mVfor 400 ms). a, control; b, during 
exposure to 100 uM UTP; c, after exposure to 
10 uM Oxo-M. eP2YP 2 R, cells expressing endog- 
enous P2Y 2 Rs. (B) Summary of Fura-4F changes 
in response to UTP (n = 18) and Oxo-M (n = 15). 
(C) Ratio of peak Fura-4F changes elicited by 
UTP and Oxo-M in nine cells. (D) Duration 
from reaching half-maximal amplitude calcium 
response to the time point of falling below 10% 
of the amplitude (n as in B). (E) Summary of 
the time to half-maximal calcium response ( n as 
in B). (F) KCNQ2/3 current inhibition in per- 
forated-patch configuration by UTP (n = 9) and 
Oxo-M (n = 7). 
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Protein lysates from tsA-201 cells probed with a pri- 
mary antibody against P2Y 2 R revealed two bands, one at 
a molecular mass of ^50 kD and the second at ^100 kD 
(Fig. 2 A) . We assume that they represent receptor mono- 
mers (predicted protein M r ^41 kD) and dimers. Lysates 
from populations of cells transfected with human P2Y 2 R 
showed greatly increased P2Y 2 R expression (both bands) 
and had to be diluted 100 times to produce bands of 
similar intensity to endogenous P2Y 2 R. We conclude that 
on average, transient overexpression increases P2Y 2 R 
by 100-fold. Presumably it is increased even more than 
the 100-fold mean in some individual well-transfected 
cells (see Falkenburger et al., 2010a). Transfecting siRNA 
against P2Y 2 R decreased the intensity of the endoge- 
nous receptor bands by 90% (Fig. 2 A). The siRNA was 
also effective at reducing the augmented expression of 
transfected P2Y 2 R RNA (Fig. 2 A) . Likewise, compared 
with control (Fig. 2 B), the P2Y 2 R siRNA decreased the 
peak amplitude of the UTP-evoked calcium rise by 90% 
in cells expressing only endogenous receptors, and 
it decreased the peak by 65% in P2Y 2 -transfected cells 
(Fig. 2, C-E). 100 pM suramin, a P2Y 2 R antagonist, re- 
duced the UTP-evoked calcium response by 68% (Fig. 2 F; 
n = 5). These separate lines of evidence indicate that 
the endogenous receptor activated by UTP is P2Y 2 . Here- 
after, we refer to endogenous P2Y 2 Rs as (low-density) 
"eP2Y 2 R," and the combination of overexpressed and en- 
dogenous P2Y 2 Rs as (high-density) "oP2Y 2 R". 



High-density P2Y 2 R can evoke PIP 2 depletion 
and KCNQ2/3 current inhibition 

To understand the apparent differences between P2Y 2 R 
and MiR actions, we next compared their effects on 
plasma membrane PIP 2 . The PIP 2 levels were monitored 
as FRET between fluorescently tagged PIP 2 -binding 
pleckstrin homology (PH) domains from PLCSj (van der 
Wal et al., 2001). Binding to PIP 2 at the plasma mem- 
brane brings CFP-PH and YFP-PH close enough together 
for FRET to occur. When PIP 2 is depleted, the PH domains 
translocate to the cytosol, and as the average distance 
between them increases, FRET decreases (Fig. 3 A) . The 
first experiments compared high-density MiR with low- 
density eP2Y 2 R. Similar to our previous studies (Jensen 
et al., 2009; Falkenburger et al., 2010a) , FRETr from the 
PH domain reporter of PIP 2 fell to ^50% of the initial 
value with the application of 10 pM Oxo-M (Fig. 3, B-D) , 
indicating a strong depletion of PIP 2 . In contrast, FRETr 
did not decrease significandy upon the application of 
100 pM UTP (Fig. 3, B and D). Thus, we observed simi- 
lar calcium responses to eP2Y 2 R and MjR activation but 
saw PIP 2 depletion and KCNQ2/3 inhibition only with 
MiR activation. We infer that the significant calcium rise 
observed with UTP neither suffices to induce net PIP 2 
depletion (by, for example, activating PLCS) nor to in- 
hibit KCNQ2/3 current by itself. 

Because of the difference in receptor densities be- 
tween MjR and eP2Y 2 R, we hypothesized that a higher 
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Figure 2. P2Y 2 is the endogenous 
receptor activated by UTP. (A; top) 
Western blot of lysates from untrans- 
fected cells, cells transfected with P2Y 2 , 
and cells transfected with siRNA against 
P2Y 2 . In the last two lanes, the lysate is 
diluted 100-fold. (Bottom) Summary 
of band densities normalized to the 
endogenous levels of P2Y 2 R from n = 4 
experiments as in A. (B) Representative 
calcium responses to 100 pM UTP and 
10 pM Oxo-M in a Fura-4F-AM-loaded 
cell without a patch pipette. (C) Rep- 
resentative calcium responses of a cell 
transfected with siRNA directed against 
the eP2Y 2 R. (D) Representative cal- 
cium responses of a cell cotransfected 
with hP2Y 2 R and siRNA against P2Y 2 R 
(E) Summary of UTP-induced calcium 
responses after siRNA transfection ( n = 6 
for control, n = 8 for siRNA, and n = 5 
for overexpression and siRNA) . (F) Cal- 
cium responses in a cell exposed to the 
P2Y 2 R antagonist suramin (100 pM). 
Representative of n = 1 cells. 
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density of activated receptors is needed to deplete PIP 2 
and to inhibit KCNQ2/3 currents than to induce cal- 
cium release. To test this hypothesis, we modified re- 
ceptor density and receptor occupancy independendy. 
First, we increased the total density of P2Y 2 R by over- 
expression, which, based on Western blot analysis, in- 
creased the density at least 100-fold (Fig. 2 A). In 
such cells with high-density P2Y 2 R, 100 pM UTP did 
reduce PH domain FRETr (Fig. 3, C and D) and in- 
hibit KCNQ2/3 current (Fig. 3, E and F) to a similar 
extent as Oxo-M acting through high-density MjR. Sec- 
ond, we lowered the concentration of Oxo-M to re- 
duce MjR occupancy. Even a very low concentration 
of Oxo-M (1 nM) elicited a Fura-4F signal of similar 
size and duration as 10 pM Oxo-M (Fig. 4, A and B), 
although the time to half-maximum was longer with 
low Oxo-M concentrations. Previous work showed no 
reduction in PH-domain FRETr and no inhibition of 
KCNQ2/3 with 1 nM Oxo-M (Jensen et al., 2009) . Thus, 
with high-density MiRs, more agonist is required to 
inhibit KCNQ2/3 and to deplete PIP 2 than to release 
calcium. In pharmacological terms, the 50% effective 
concentration for calcium release is much lower than 
that for KCNQ2/3 inhibition; i.e., the number of 



spare receptors, the receptor reserve, is large for cal- 
cium release and smaller for KCNQ2/3 inhibition and 
PIP 2 depletion. 

These findings demonstrate that agonist and recep- 
tor requirements can differ between different outputs 
of G q -coupled signaling and suggest the hypothesis 
that the differences we see between P2Y 2 R and MiR 
actions might be explained entirely by a difference in 
receptor density. Our companion paper (Falkenburger 
et al., 2013) describes a mathematical, kinetic model 
based on this hypothesis that reproduces most of the 
features we have seen. Figs. 4 (C and D) and 5 (A and B) 
in this paper show simulations from that model. With 
appropriate choices of parameters, an increase in recep- 
tor density from 1 to 500 per pm 2 converts a response 
with little depletion of PIP 2 or decrease of current, yet still 
a strong Ca 2+ rise and Fura-4F response (Fig. 5, A and B, 
"UTP") , into a response with the typical full Fura-4F sig- 
nal, full inhibition of KCNQ current (Fig. 5 A "Oxo-M") , 
strong Ca 2+ increase, and full depletion of PIP 2 (Fig. 5 B, 
"Oxo-M"). The same calculations show that very low con- 
centrations of agonist still evoke a full-amplitude Fura-4F 
response and a strong Ca 2+ transient when receptor den- 
sity is high (Fig. 4, C and D). 
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Figure 3. Activation of endogenous P2Y 2 R does not pro- 
duce a net reduction in PIP 2 . (A) Schematic representation 
of PIP 2 measurements by PH-domain FRET. Binding of 
PH-CFP and PH-YFP to PIP 2 at the plasma membrane brings 
CFP and YFP close enough together for FRET. After PIP 2 
hydrolysis by PLC, PH-CFP and PH-YFP move apart into the 
cytosol. (B) FRETr time course from a single cell expressing 
PH-CFP and PH-YFP during exposure to 100 uM UTP and 
10 uM Oxo-M. (C) FRETr time course from a cell express- 
ing hP2Y 2 R MjR, CFP-PH, and YFP-PH. (D) Summary of 
percent FRETr change in experiments like those in B and C 
(n = 10 with endogenous P2Y 2 , n = 12 with high-density 
transfected P2Y 2 , and n = 13 with high-density transfected 
MjRs). (E) Representative KCNQ2/3 tail-current changes 
with UTP in a cell transfected with KCNQ2, KCNQ3, and 
P2Y 2 R. (F) Summary of KCNQ inhibition in experiments 
like those shown in E and Fig. 1 A (n = 1 1 for endogenous 
P2Y 2 and n = 5 for high-density P2Y 2 ). 
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Figure 4. Similar calcium responses are evoked by low and high Oxo-M concentrations. (A) Representative calcium responses of a cell 
transfected with IVT^Rs, loaded with Fura-4F-AM, and treated with 1 nM, 0.1 pM, or 10 pM Oxo-M. (B) Summary of Fura-4F amplitude, 
time to half-maximum, and duration for experiments as in A (n = 16) . (C) Simulations from our kinetic model reproducing the observa- 
tions in A. Quotation marks on Oxo-M are a reminder that this is a computer simulation rather than experiment. Assumed [Fura-4F] was 
1 pM. (D) Simulation showing corresponding rises in cytosolic calcium concentration during 1 nM, 0.1 pM, and 10 pM Oxo-M. Fura-4F 
concentration was set at 0 pM. 



Receptors stimulate IP 3 production 

To identify at which stage of the G q signaling cascade 
the receptor reserve arises, we measured IP 3 production 
in living tsA-201 cells using two FRET reporters, IRIST 
(Matsu-ura et al., 2006) and LIBRAvIII (Tanimura et al., 
2009) . Both consist of CFP and YFP linked by the ligand- 
binding domain of the rat IP 3 R type III. They show a 
reduction in FRET between CFP and YFP upon binding 
of IP 3 . IRIST is a cytosolic protein, whereas LIBRAvIII is 
targeted to the plasma membrane by a palmitoylation 
sequence from GAP43 (Fig. S3). The reported ligand con- 
centrations for half-maximal FRET change are similar: 



550 nM IP 3 for IRIS-1 and 490 nM for LIBRAvIII. These 
probes are insensitive to physiological changes of cal- 
cium. In our hands, IRIS-1 signals were difficult to inter- 
pret and were not graded with Oxo-M concentration. 
On average, the IRIS-1 response is an ^15% FRET de- 
crease with invariant temporal characteristics in response 
to 1 nM, 100 nM, and 10 uM Oxo-M (Fig. S4, A and B, 
and not depicted). Further, with 10 uM Oxo-M, the 
delay before onset, the onset and recovery time constants, 
and the IRIS-1 response duration were longer than ex- 
pected from our calcium measurements or from re- 
ported biochemical assays of IP 3 (Willars et al., 1998), 
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Figure 5. Receptor density can 
account for differences in calcium 
signaling and PIP 2 depletion. 
Simulations from our kinetic 
model. Low receptor densities 
(1 per pm 2 ) mimic eP2Y 2 , and high 
receptor densities (500 per pm 2 ) 
mimic high-density M^. (A) Sim- 
ulated time courses of Fura-4F 
(top) and KCNQ2/3 current (bot- 
tom) responses to UTP acting on 
low receptor densities (red line) 
and Oxo-M acting on high re- 
ceptor densities (black line). 
[Fura-4F] was 1 pM. (B) Simul- 
ated time courses of calcium 
(top) and PIP 2 (bottom) concen- 
tration changes with low and high 
receptor densities. Fura-4F con- 
cenuation was set at 0 pM. 
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and much longer than those observed with LIBRAvIII. 
Possible explanations for such unexpected kinetics in- 
clude reporting of IP3 in cellular regions far from the 
plasma membrane. 

We therefore turned to membrane-localized LIBRAvIII. 
It responded to 10 pM Oxo-M with a small, reproduc- 
ible FRETr decrease (Fig. 6 A) . The IP3 selectivity was 
assessed by overexpressing the enzyme IP3 5-phospha- 
tase, which can deplete IP3 as fast as it is made in these 
cells (Horowitz et al., 2005). Overexpression of this en- 
zyme eliminated the responses of LIBRAvIII (Fig. 6 A), 
confirming the FRET probe's selectivity for IP 3 . We 
calibrated the probe by dialyzing different concentra- 
tions of IP3 into the cell via a patch pipette (Fig. 6 B, 
circles). Time-dependent FRET changes began in <10 s 
after breakthrough. No detectable change in FRETr oc- 
curred with dialysis of 1 pM IP3, whereas the response 
nearly saturated with dialysis of 10 pM IP3. Thus, the 
LIBRAvIII probe has a fairly narrow dynamic range. To 



understand the dialysis protocol better, the entry of 
IP3 was simulated by adding to our larger cell model 
a first-order exchange from the pipette through a rate- 
limiting orifice. This model reproduced the calibra- 
tion experiment reasonably well (lines in Fig. 6 B), 
assuming a dissociation constant of 500 nM for the 
IP3-LIBRA complex as reported by others (Tanimura 
etal., 2009). The model indicates that the endogenous 
steady-state IP3 5-phosphatase activity, which allows a 
mean lifetime for IP 3 of only 12.5 s, severely reduces 
the cellular IP 3 during dialysis. Separate control experi- 
ments (not depicted) with fluorescent dyes in the pi- 
pette typically showed dialysis of dye into the cytosol 
with an exponential time constant of 50-70 s. With sim- 
ilar pipette exchange rates for IP 3 added to the model 
and retaining an active endogenous IP3 phosphatase, 
the steady-state cytosolic IP3 concentration reaches only 
roughly 20% of that in the pipette, and the time con- 
stant of the exponential change is only 20% of that with 
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Figure 6. LIBRAvIII reports changes in 
IPs concentration before a rise in calcium. 
(A) Averaged time courses of normalized 
FRET from IP 3 probe LIBRAvIII in response to 
10 pM Oxo-M. Cells were transfected with 
IVLR, LIBRAvIII, and with (n = 6) or without 
(n= 7) IP 3 5-phosphatase. Points were mildly 
filtered (see Materials and methods). (B) Cal- 
ibration of LIBRAvIII. Cells were transfected 
with LIBRAvIII and patched with pipettes 
containing different concentrations of IP3 
(1-100 pM). Plasma membrane was ruptured 
at t = 0 s, leading to dialysis of IP 3 into the cyto- 
sol of the cell (lines with markers) . A response 
to 1 0 pM Oxo-M is superimposed for compari- 
son. Smooth curves without markers are model 
simulations of IP 3 entry assuming 6 pM of 
cytosolic LIBRAvIII, a dissociation constant 
for LIBRAvIII of 0.5 pM, and a 12.5-s lifetime 
for IP 3 in the cytosol (see Falkenburger et al., 
2013 for details). (C) Representative FRET re- 
cording of a cell expressing LIBRAvIII and 
MjR exposed to a series of increasing Oxo-M 
concentrations. Points were mildly filtered 
(see Materials and methods). (D) Summary of 
LIBRAvIII responses to 0.1 and 10 pM Oxo-M 
(n = 7). (E and F) Simulations showing the per- 
centage of IP3 bound to LIBRAvIII (E) and 
the concentration of IP 3 (F) after activation of 
MjR (density, 500 per pm 2 ) with 1 nM, 100 nM, 
or 10 pm Oxo-M. [LIBRAvIII] was 6 pM for E 
and 0 pM for F. (G and H) Superimposed nor- 
malized experimental Fura-4F and LIBRAvIII 
responses to 10 pM Oxo-M (n = 7). G shows 
an expanded version of the dashed rectangle 
in H. Note that the scales for LIBRAvIII in 
G and H are inverted. 
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no cytosolic phosphatase (see Fig. S5 A) . In addition 
LIBRAvIII significantly buffers the IP 3 (see Fig. S5 B). 

LIBRAvIII responded to G q PCR activation in a dose- 
dependent manner. No response was detected with 
1 nM Oxo-M, whereas with 100 nM or 10 pM Oxo-M the 
amplitudes were not statistically different, which is con- 
sistent with the narrow dynamic range observed during 
IP3 dialysis. However, the time courses differed (Fig. 6, 
C and D) . Comparing responses to 100 nM versus 10 pM 
Oxo-M, the delay before onset was longer (18 ± 2 s vs. 
9 ± 6 s) , the onset slower (T on 7 ± 1 s vs. 3 ± 1 s) , and the 
duration shorter (40 ± 3 s vs. 72 ± 10 s; Fig. 6 D). The 
slower response and faster recovery are expected from 
the weaker activation of PLC by 100 nM Oxo-M (Jensen 
etal., 2009). 

What is the range of IP 3 changes? When the Oxo-M 
responses of LIBRAvIII were compared with those in 



the IP 3 calibration experiments, we found that 10 pM 
Oxo-M elicited responses similar in amplitude and 
kinetics to dialysis with 10 pM IP 3 in the pipette. The 
IP 3 -binding curve of LIBRAvIII (Fig. S5 C) shows that 
LIBRAvIII does not discriminate well between IP 3 con- 
centrations above 2 pM. We conclude that peak IP 3 is 
above 2 pM. The resting levels of IP3 are below the 
level reached with 2 pM IP 3 in the pipette (<300 nM 
in the cell), which evoked a minimal FRETr decrease. 
Fig. 6 (E and F) shows that the dependence of the 
LIBRAvIII response amplitude, rise time, and dura- 
tion on Oxo-M concentration was captured qualita- 
tively in the kinetic model described in our companion 
paper (Falkenburger et al., 2013). In this model, free 
IP 3 rises to 10 pM in the absence of LIBRAvIII. This 
limit would result from cleavage of nearly all of the cel- 
lular PIP 2 . 
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Figure 7. A little IP 3 suffices to initiate robust 
calcium signals. (A) Time course of calcium 
responses (Fura-4F) to UTP application. Cells 
were transfected with plasma membrane- 
targeted FRB and the PIP 2 5-phosphatase 
CFP-FKBP-Inp54p for recruitment of the 
5-phosphatase to the plasma membrane by 
5 uM rapamycin (bar). Representative of n = 5 
cells. (B) Simulations from our kinetic model 
of the observations in A. A second calcula- 
tion starts at 500 s in which a supplemental 
5-phosphatase activity is turned on. Note that 
as PIP 2 gradually falls to a new steady state, 
less resting IP 3 is made by the basal PLC activ- 
ity (a function of basal nucleotide exchange 
at Get), and basal calcium release declines, 
slowly decreasing resting Fura-4F saturation. 
(C) Simultaneous recording of KCNQ2/3 
tail current (bottom) and cytosolic calcium 
by Fura-4F ratio (top). Cells were transfected 
with MiR, KCNQ2/3 channels, and the VSP. 
During baseline, KCNQ2/3 channels were 
activated maximally by depolarization from 
-60 to +40 mV for 200 ms every 600 ms. Dur- 
ing the shaded rectangle marked "VSP," de- 
polarization to +120 mV activated KCNQ2/3 
channels as well as VSP, which depleted PIP 2 
by ~80% within seconds. Then, 10 uM Oxo-M 
was applied. The resulting calcium response 
was compared with calcium responses elicited 
before VSP-induced PIP 2 depletion. (D) Sim- 
ulations showing Fura-4F saturation with cal- 
cium (top) and KCNQ2/3 current (bottom) 
in response to 10 pM Oxo-M at t = 0 s (black) 
and with activation of VSP from —50 s until 
150 s (red). The extent of VSP activation was 
adjusted to have ^90% KCNQ2/3 current in- 
hibition (VSP_max = 0.3 s _1 ). Note the reduc- 
tion of the duration but not the amplitude of 
the Fura-4F signal by VSP. (E and F) Summary 
of experiments as in C. The recorded cells are 
grouped by extent of KCNQ2/3 inhibition 
(percent PIP 2 depletion) for summaries. 
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A low threshold for IP 3 -dependent calcium release 

We now turn to Ca 2+ release from the ER. To learn 
the kinetic relation of IP3 production to calcium re- 
lease, we overlaid normalized experimental Fura-4F and 
LIBRAvIII responses for 10 pM Oxo-M (Fig. 6, G and H) . 
LIBRAvIII reports a rise in IP 3 beginning already at 
2 s, whereas the sudden upstroke of the Fura-4F sig- 
nal occurs only between 4 and 6 s (Fig. 6 G). Hence, 
LIBRAvIII is a fast reporter of IP3, and IP3 production 
starts early. However, once initiated, the resulting calcium 
response rises more rapidly and peaks sooner than the 
IP3 signal (Fig. 6 G). The Ca 2+ rise, in particular its ste- 
reotypical amplitude and the concentration-dependent 
time to half-maximum (Fig. 4 B), is nearly like a regen- 
erative response, rising abrupdy with positive feedback 
when graded IP3 accumulation reaches a certain level. 
The decay time courses for calcium and IP3 are quite dif- 
ferent (Fig. 6 H), reflecting their mechanistically sepa- 
rate clearance mechanisms. 

Despite the strong Fura-4F signal produced by 1 nM 
Oxo-M (Fig. 4 A) , the amount of IP3 made at that low 
agonist concentration is too small to be reported by 
LIBRAvIII (Fig. 6 C) . A minimal production of IP3 suf- 
fices to initiate the calcium transient. This conclusion 
was reinforced by experiments showing little change 
of calcium responses upon reducing the available PIP2 
with coexpressed PIP2 5-phosphatases. Such prior PIP 2 
depletion ought to reduce the amount of IP3 that can 
be produced by a given stimulus. First, we used rapamy- 
cin-induced dimerization of an FKBP domain with an 
FRB domain to "activate" the PIP 2 5-phosphatase Inp54p. 
When FKBP-tagged Inp54p is expressed together with 
plasma membrane-targeted FRB, the addition of rapa- 
mycin leads to translocation of the phosphatase to the 
plasma membrane and strong PIP 2 depletion (Suh et al., 



2006). However, this manipulation did not reduce the 
amplitude of subsequent calcium responses to UTP 
(Fig. 7 A). The model also shows a resistance of the 
calcium response to 90% depletion of PIP 2 (Fig. 7 B). 
A second approach used a VSP to deplete PIP 2 and 
gave the same result (Fig. 7 C). Depolarization activates 
VSP reversibly and depletes plasma membrane PIP 2 
(Murata and Okamura, 2007; Falkenburger et al., 2010b) . 
PIP 2 levels were continuously monitored by recording 
KCNQ2/3 currents. To record KCNQ2/3 currents while 
keeping PIP 2 levels low over tens of seconds in VSP- 
transfected cells, we stepped cyclically from a holding po- 
tential of -60 to +120 mVfor 400 ms, every 800 ms, and 
measured KCNQ2/3 tail currents at —60 mV. On aver- 
age, the intermittent depolarization to +120 mV was suf- 
ficient to reduce KCNQ2/3 current by ^80%. The time 
course of PIP 2 reduction as reported by PH-domain FRET 
was saw-toothed with an overall relaxation time con- 
stant of 10 s (Fig. S6). Despite this accumulating PIP 2 
depletion, 10 pM Oxo-M could still evoke a full-ampli- 
tude calcium response (Fig. 7 C) . Because the KCNQ2/3 
current reduction achieved by the intermittent depolar- 
ization to +120 mV varied with VSP expression levels, we 
grouped the recorded cells by the extent of KCNQ2/3 
inhibition for quantification (Fig. 7, E and F) . Simula- 
tions from our kinetic model show that when VSP activ- 
ity was adjusted to give a 90% reduction in KCNQ2/3 
current, 10 pM Oxo-M generated a calcium response of 
similar amplitude but reduced duration (Fig. 7 D) , simi- 
lar to our experimental data (Fig. 7, E and F) . In sum- 
mary, activation of PLC can generate sufficient IP3 for a 
calcium response even when PIP 2 levels are only a frac- 
tion of normal. It should be emphasized that the two 
5-phosphatase treatments we used do lower PIP 2 consider- 
ably but never to zero because PIP 2 synthesis continues. 
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Figure 8. PIP 2 synthesis is required 
for maintained IP 3 production. To 
determine whether PIP 2 is synthesized 
during G q activation, we inhibited PI 
4-kinases by 30 uM wortmannin dur- 
ing a prolonged application of Oxo-M 
(200 s at 10 uM) and measured IP 3 
production by LIBRAvIII. Cells were 
exposed to wortmannin for 5 min be- 
fore start of recording. (A and B) Time 
course of normalized LIBRAvIII re- 
sponses (control, n = 8 cells; wortman- 
nin, n = 5). (C and D) Simulations 
from model. In C, the PI 4-kinase re- 
action was accelerated during agonist 
application with the "standard" time 
course shown in Fig. S9 A. In D, the PI 
4-kinase was either left unchanged at its 
resting rate or turned off fully so that 
only the PIP 2 already present at rest is 
available to generate IP3. To achieve a 
stable baseline with PI 4-kinase off, the 
PI 4-phosphatase reaction was turned 
off as well in this case. 
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Ongoing PIP 2 synthesis is required for continual 
IP 3 production 

Is there significant synthesis of PIP 2 during agonist 
application? Throughout a 200-s application of 10 uM 
Oxo-M, IP 3 remained elevated, with some possible sag 
after an initial peak (Fig. 8 A) . Given that (a) the life- 
time of IP3 is short (Fig. 6) and (b) PIP 2 is depleted by 
90% during 10 uM Oxo-M (Horowitz et al., 2005), this 
continued production of IP3 implies that new precursor 
PIP 2 is being formed as PLC is breaking it down. The 
rate-limiting step in PIP 2 synthesis is the PI 4-kinase 
activity; the PIP 5-kinase activity that makes PIP 2 from 
PI(4)P is 20 times faster (Falkenburger et al., 2010b) . The 
slow step is probably accelerated during agonist appli- 
cation (Xu et al., 2003; Falkenburger et al., 2010b). 
To determine the importance of de novo PIP 2 synthesis 
during Oxo-M stimulation, we inhibited the PI 4-kinase 
using 30 |iM wortmannin. After wortmannin, the initial 
amplitude of the LIBRAvIII FRETr response to Oxo-M 
was normal, but then it relaxed back to baseline while 
Oxo-M was still present (Fig. 8 B): IP 3 production had 
stopped. Fig. 8 (C and D) shows that a model with ac- 
celerated PIP 2 synthesis during Oxo-M application can 
explain maintained IP3 elevation throughout Oxo-M ap- 
plication, and that continued IP3 production fades if PI 
4-kinase is not accelerated or is turned off fully. 

DISCUSSION 

To provide a more complete description of signaling by 
two different G q -coupled receptors, we have performed 
quantitative kinetic measurements of PIP 2 , IP3, calcium, 
and KCNQ2/3 channel current in response to agonist 
application. Our measurements in tsA-201 cells show 
that (a) endogenous P2Y 2 Rs activate the PLC pathway 
and signal to KCNQ channels much less effectively than 
MiRs expressed at high density; (b) experimentally, these 
differences are largely removed when P2Y 2 Rs are ex- 
pressed at high density; (c) very little IP3 needs to be gen- 
erated to induce a robust calcium release; (d) less receptor 
stimulation is needed to elicit a Ca 2+ elevation than to 
deplete PIP 2 or to inhibit KCNQ channels; and (e) accel- 
erated PIP 2 synthesis probably boosts IP 3 production 
during long receptor activation. These measurements fur- 
ther guided and are explained by an expanded kinetic 
model detailed in our companion paper (Falkenburger 
et al., 2013). 

Differences between P2Y 2 Rs and MiRs 
We first revisit the contrasting actions of endogenous 
P2Y 2 Rs and high-density MiRs. A saturating concentra- 
tion of UTP acting on eP2Y 2 Rs leads to calcium release 
without inhibition of KCNQ2/3 channels, whereas a sat- 
urating concentration of Oxo-M acting on high-density 
MjRs does both. Nevertheless, we conclude that there is 
no qualitative difference between P2Y 2 R and MiR actions. 



Our measurements show that endogenous P2Y 2 Rs are 
at least 100-fold less abundant than overexpressed 
P2Y 2 Rs, and our modeling shows that a 100-500-fold dif- 
ference in receptor density between endogenous P2Y 2 Rs 
and overexpressed MiRs is sufficient to explain all ob- 
servations. In support, decreasing the number of acti- 
vated MiRs by using an extremely low concentration 
of agonist ( 1 nM Oxo-M) mimics the responses to endog- 
enous P2Y 2 R activation by a saturating concentration of 
UTP. Conversely, increasing the number of P2Y 2 Rs 
by overexpression makes the effects of UTP mimic those 
of Oxo-M. The density of endogenous purinergic re- 
ceptors that we estimate is in a range expected from 
our previous estimate of endogenous G proteins and a 
typical ratio of G proteins to receptors (Falkenburger 
et al., 2010a). We therefore attribute the observed dif- 
ferences between UTP and Oxo-M primarily to a differ- 
ence in receptor numbers. Such observations illustrate 
how quantitative differences in the abundance of a sig- 
naling component can result in apparent qualitative dif- 
ferences in responses. 

Comparison with neurons 

In sympathetic neurons, activation of MjRs leads to PIP 2 
depletion (Zaika et al., 2011), and activation of B 2 re- 
ceptors leads to calcium release (Delmas and Brown, 
2002). Based on the findings described here, the first ob- 
servation might be explained by a high density of MjRs 
and a low density of B 2 receptors in these neurons. Of 
course, receptors can be arranged in more interest- 
ing ways than just being high or low density, particu- 
larly in differentiated cells with specialized membrane 
compartments (e.g., spines, cilia, microvilli, or the apical 
membranes of epithelia) . What counts for downstream 
effects is not the overall density of receptors but the 
local density "seen" by PLC molecules. Our data in tsA201 
cells suggest that we do not need many PLC molecules 
(10 per um 2 ) to fully deplete PIP 2 with the time course and 
concentration dependence observed experimentally. 
The interaction of signaling molecules will be strongly 
affected by inhomogenous densities resulting from mem- 
brane compartmentalization and scaffolding molecules. 
Nonetheless, we believe that the estimates for affinities 
and binding kinetics between MjR, G q , and PLC derived 
from our previous experiments and modeling (Jensen 
et al., 2009; Falkenburger et al., 2010a, b) are accurate 
and similarly apply to sympathetic neurons. Consequently, 
we have to assume at least 10 pm molecules of PLC 
to deplete PIP 2 in sympathetic neurons and a density of 
MjR that is high, at least in the vicinity of PLC. 

It is not unexpected that low-density B 2 receptors 
can lead to calcium release, but it is hard to explain why 
MjRs do not. Others have suggested that B 2 receptors 
are close to IP 3 Rs and MiRs are not, and that fast IP 3 deg- 
radation prevents its diffusion from MiR-activated PLC 
to IP 3 Rs (Delmas et al., 2004). Expanding measurements 
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and modeling into these spatial aspects will be an im- 
portant next step to further understand how signaling 
specificity of GqPCR works. 

Calcium signaling needs only minimal stimulation of PLC 

Different signaling end points require different inten- 
sities of stimulation of receptors and PLC. For exam- 
ple, MjR agonist binding has a midpoint of 4-10 pM for 
Oxo-M (Jensen et al., 2009), yet with high-density MiR, 
closure of KCNQ2/3 channels and FRET interactions 
between Ga q and PLC and those between M4R and G 
all have midpoints at 120-330 nM Oxo-M (Jensen et al., 
2009) . The release of stored calcium by 1P3 seems to be 
engaged by the smallest stimulation of receptors and 
PLC in our tsA cells: (a) Exposure to 1 nM Oxo-M acti- 
vates PLC sufficiendy to evoke a calcium rise, even though 
it produces too little IP3 to be detected by the LIBRAvIII 
IP3 reporter, (b) Calcium is released by activation of 
endogenous UTP receptors without a measurable net 
depletion of plasma membrane PIP 2 . (c) A significant 
calcium release can be elicited even after depleting much 
of the normal PIP 2 using exogenous PIP 2 5-phosphatases, 
whether by rapamycin-induced protein dimerization 
or by VSP combined with membrane depolarization. 
These manipulations reduced PIP 2 , and probably IP3 pro- 
duction, by up to 90% as judged from KCNQ2/3 current 
inhibition. In drug discovery with overexpressed G q PCRs, 
calcium measurements are often used as an assay for 
receptor activation. Such assays should give extremely 
sensitive responses. The release takes on a somewhat 
all-or-none character from regenerative positive feed- 
back and occurs at the very low end of receptor occu- 
pancy (see below) rather than being broadly graded with 
receptor activation. 

Several factors conspire to allow a small stimulus to 
make a strong calcium response and for the calcium re- 
sponse to be less graded as a function of agonist con- 
centration. We use ideas from published experiments 
corroborated by our modeling (see Falkenburger et al., 
2013) as an explanation. As is discussed in our compan- 
ion paper, there are many observations of the IP3 and 
calcium dependence of IP3R channel opening. They 
all agree that both IP 3 and calcium show cooperativity, 
with channel activation following power law depen- 
dences (second-fourth power) of their concentrations 
that steepen the activation curve, as originally described 
by Bezprozvanny et al. (1991) and reviewed by Foskett 
et al. (2007). Calcium also activates PLC, thus accelerat- 
ing IP3 production. In our model, IP3 turns on calcium 
release with as litde as 1% full IP3 production (which 
can rise as high as ^12 pM IP 3 for the strongest acti- 
vation). Half-maximal activation of channel opening 
by 100-300 nM Ca 2+ is augmented by local microdo- 
mains of calcium rise around active groups of receptors. 
Finally, in real cells, but not part of our model, ryanodine 
receptors on the ER can also add to the Ca 2+ -induced 



Ca 2+ release. There are conflicting reports in the literature 
concerning ryanodine receptors, but our cells responded 
robusdy with a calcium transient upon application of 
caffeine (Fig. S7), suggesting that tsA-201 cells do express 
ryanodine receptors and that a more complete model 
should include them. Hence, although the time course 
and duration of calcium signals is somewhat graded 
with stimulus intensity, the peak amplitude rises fairly 
abrupdy, and the range over which the response ampli- 
tude is graded is compressed. 

There are many published partial models of the 
IP3R. The rationale for design and parameters of our 
model are detailed in our companion paper. We now 
consider how our model explains calcium transients in 
terms that are instructive but quite likely would differ 
in detail had we chosen another model of the IP3R: 
The calcium response even to short applications of 
10 pM Oxo-M is protracted in a late "hump" or plateau 
(Figs. 2, C and D, and 4 A) . During longer Oxo-M ap- 
plications (200 s), the response could be reduced and 
shortened by switching from 2 to 0 mM of extracellu- 
lar calcium (see Fig. 9 in Horowitz et al., 2005), a stan- 
dard test for store-operated calcium entry (SOCE). Our 
model does not include SOCE. Nevertheless, it does 
produce a plateau that prolongs the calcium response 
(Figs. 5 A and 8, B and C) . This plateau results from 
unanticipated delayed and nonlinear effects of cal- 
cium binding to the IP3R. In the model, >99.7% of the 
IPsRs become inactivated within the first second as 
a consequence of rising calcium (Fig. S8 A, 1-h 3 ). Cal- 
cium is then gradually cleared through the SERCA 
pump, and as calcium falls again, a small percentage 
(0.1%) of receptors recover from inactivation by t = 20 s 
(while IP3 is still high) , producing a second phase of 
calcium flux in IP 3 Rs. The resulting calcium plateau 
represents a dynamic steady state of calcium release 
and calcium clearance, and its height in our model is 
determined by the maximum transport capacity of the 
SERCA pump (vP; Fig. S8 B). The subsequent, final 
decrease in cytosolic calcium is caused by the fall in 
occupancy of IP 3 -binding sites on the IP3R that occurs 
with declining IP 3 levels. The duration of the calcium 
response is a function of the apparent affinity of the 
IP3R for IP3 and was determined primarily by the dura- 
tion of the IP3 elevation and by Ki P3 (Fig. S8 C) . We con- 
clude that a plateau in the calcium response can result 
from intrinsic properties of the IP3R and is not neces- 
sarily a sign of SOCE. 

In pharmacological terms, there are spare receptors, 
a significant receptor reserve, for calcium signaling. Thus, 
according to our model, 1 nM Oxo-M would occupy 
only 0.06% of the overexpressed MiRs and activate 
only 0.3% of the endogenous PLC. As a consequence, 
0.1% of the PIP 2 is hydrolyzed per second. In the ab- 
sence of LIBRAvIII, cytosolic IP3 would rise to 200 nM 
after 20 s. This amount is well above the half-activating 
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concentration for IP3R calcium release yet below the 
500-nM half-saturation concentration for detection by 
the lower-affinity LIBRAvIII reporter. 

PIP 2 synthesis must be accelerated during 
GqPCR activation 

A couple of seconds after agonist addition, LIBRAvIII 
starts reporting an increase in IP 3 , and the IP3 elevation 
lasts throughout a 200-s agonist application despite the 
major decline of PIP 2 concentration and KCNQ cur- 
rents. The IP3 elevation requires continued PIP 2 synthe- 
sis, as it is not maintained after wortmannin treatment, 
which blocks type III PI 4-kinase. Similar persistent IP3 
production, sensitive to wortmannin, has been shown in 
experiments with radiolabeled inositol in SH-SY5Y cells 
(Nakanishi et al., 1995; Willars et al., 1998). Our stan- 
dard kinetic model has resting rates for PI 4-kinase and 
PIP 5-kinase, chosen to match the recovery of PIP 2 both 
after depletion by PLC and after depletion by VSP 
(Horowitz et al., 2005; Falkenburger et al., 2010b) . The 
5-kinase is ^20 times faster than the 4-kinase, making 
the wortmannin-sensitive PI 4-kinase the rate-limiting 
enzyme of PIP 2 synthesis. In the model, the mean life- 
time of an IP3 molecule is only 12.5 s, and we find that 
the chosen resting rates of PIP 2 synthesis are too slow 
to support the observed persistent IP3 production dur- 
ing agonist application. These observations require that 
the rate of the PI 4-kinase be elevated temporarily dur- 
ing G q PCR activation as was done in the model of Xu 
etal. (2003). 

We have not studied the biochemical mechanism 
of kinase acceleration. There are several types of PI 
4-kinase that synthesize PI(4)P. Shapiro and colleagues 
found that the inhibition of calcium currents in supe- 
rior cervical ganglion cells by G q PCR activation, which 
is attributed to PIP 2 depletion, is facilitated by dominant- 
negative neuronal calcium sensor 1 (NCST) (Camper 
et al., 2004; Zaika et al., 2007, 2011). NCS-1 binds cal- 
cium and PI 4-kinase 111(3 (Zhao et al., 2001). In PC12 
cells, overexpression of NCS-1 increased basal levels of 
PI(4)P and PI(4,5)P 2 and consequently G q PCR-induced 
IP3 production (Koizumi et al., 2002), suggesting a reg- 
ulation of constitutive PIP 2 synthesis. The inhibition of 
receptor-induced responses by wortmannin (Nakanishi 
et al., 1995) suggests a type III PI 4-kinase, and indeed 
PIK93, a specific inhibitor of PI 4-kinase III|3, facilitated 
inhibition of calcium currents in superior cervical gan- 
glion cells (Zaika et al., 2011). Yet, PIK93 did not affect 
IP3 production by angiotensin or recovery of PI(4)P 
and PI(4,5)P 2 in HEK293 cells. Instead, PI(4,5)P 2 levels 
were most strongly depressed during angiotensin ap- 
plication in cells expressing siRNA for PI 4-kinase Ilia 
(Balla et al., 2008). Earlier studies had suggested that 
neither subtype is localized to the plasma membrane 
(Balla, 2007) , but PI4-kinase Ilia was recently found 
to be targeted to the plasma membrane by a complex 



containing several additional proteins (Nakatsu et al., 
2012). As an alternative to PI 4-kinase acceleration, PI (4) 
P supply might be accelerated during G q PCR activation 
by increased delivery of PI(4)P to the plasma membrane, 
for instance, by vesicle trafficking or by the sterol transfer 
protein Osh4p (de Saint-Jean etal., 2011) from the Golgi. 

There are also several types of PIP 5-kinases — la, 1(3, 
and I7 — that synthesize PI(4,5)P 2 . All are constitutively 
inhibited by phosphorylation. They are activated by 
phosphatidic acid, the product of DAG kinase (Kanaho 
et al., 2007) , activated by the small GTPase RhoA (Oude 
Weernink et al., 2004) , which can be stimulated by G q PCR 
(Dutt et al., 2002), and activated by dephosphorylation 
through phosphatases that can be stimulated by calcium 
elevation (Unoki et al., 2012). Inhibition of Rho facili- 
tated PIP 2 depletion during G q PCR activation in supe- 
rior cervical ganglion cells (Zaika et al., 2011). The PIP 
5-kinase subtype involved in G q PCR-induced synthesis 
appears to depend on the cell type (Tolias et al., 2000; 
Wang etal., 2004, 2008). 

Although we do not know the biochemical mecha- 
nism, our model provides constraints for the time course 
of the acceleration of PIP 2 synthesis during G q PCR acti- 
vation. Assuming a very fast onset of acceleration in the 
model (red trace in Fig. S9 A) produced a transient 
increase in PIP 2 (red trace in Fig. S9 B), which we did 
not observe experimentally when recording KCNQ2/3 
current. Therefore, onset may rise with a time constant 
of 1 s or longer after receptor activation (black traces 
in Fig. S9). However, if it rises very much slower than 
this, steady state is not reached within 20 s, which also 
contradicts our experimental findings. Second, if the 
acceleration of PIP 2 synthesis decays too quickly after 
agonist wash-off, PIP 2 levels would continue to drop, 
which we also did not observe. Therefore, acceleration 
of PIP 2 synthesis may decline with a time constant of 
10 s or longer. We tested whether the model would 
work assuming that resting kinase rates are zero. It did 
not. We found that PIP 2 recovery after Oxo-M applica- 
tion occurs too quickly (green trace in Fig. S9) , over- 
shoots, or is insufficient, unless a steady resting value 
for PI 4-kinase activity is assumed (black trace in Fig. S9) . 
Although we did not include it in the model, the accel- 
eration of PIP 2 synthesis might start faster and last lon- 
ger as the concentration of Oxo-M is raised. Indeed, 
the model does not fully reproduce the observed dif- 
ference in response duration between 0.1 and 10 uM 
Oxo-M for the DAG and IP3 reporters, and the time 
constants derived for 10 pM Oxo-M produce the men- 
tioned transients for other concentrations (see Fig. S4 B 
of our companion paper) . Because these findings suggest 
that kinase acceleration tracks G q PCR activation fairly 
closely, calcium with its transient and more autonomous 
time course is perhaps not as likely to underlie the accel- 
eration of PIP 2 synthesis as, for example, G protein a 
or (37 subunits. 
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At this point, the major unresolved question for the 
future is how to introduce acceleration of the lipid ki- 
nases into the model in a proven mechanistic manner 
that depends on other components of the model rather 
than using an empirical ad hoc time course of accelera- 
tion. Such a model should include in a natural way any 
dependence of acceleration on time and on the strength 
of receptor activation. In addition the complexities of 
calcium dynamics, clearance, buffering, and SOCE re- 
main to be included. 

Conclusion 

Our experiments and modeling have shown that appar- 
ent qualitative differences in downstream signaling from 
muscarinic and purinergic G q PCRs can be understood 
as entirely caused by differences in receptor density. 
We also find that a very low threshold for IP 3 -dependent 
calcium release and the cooperative nature of the release 
kinetics mean that large, nearly regenerative calcium 
signals are elicited by quite weak stimulation of G q PCR 
pathways. Finally, we find that receptor stimulation is ac- 
companied by acceleration of the enzymes of PIP 2 syn- 
thesis, but we do not know the underlying mechanism. 
Our companion paper continues the analysis of DAG 
production and PKC stimulation and gives further details 
of our model. 
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